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Abstract 



We describe an additional module for the Mathematica package Feyn- 
Rules that allows for an easy building of any N = 1 supersymmetric quan- 
tum field theory, directly in superspace. After the superfield content of a 
specific model has been implemented, the user can study the properties of 
the model, such as the supersymmetric transformation laws of the associ- 
ated Lagrangian, directly in Mathematica. While the model dependent 
parts of the latter, i.e., the soft supersymmetry-breaking Lagrangian and 
the superpotential, have to be provided by the user, the model independent 
pieces, such as the gauge interaction terms, are derived automatically. Using 
the strengths of the FeynRules program, it is then possible to derive all 
the Feynman rules associated to the model and implement them in all the 
Feynman diagram calculators interfaced to FeynRules in a straightforward 
way. 
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1. Introduction 

Supersymmetric (SUSY) theories are among the most popular extensions 
of the Standard Model (SM) of particle physics and widely studied both 
on the theoretical and experimental sides. Apart from relating bosons with 
fermions and unifying internal and space-time symmetries, they address a 
set of conceptual problems of the SM, such as the large hierarchy between 
the electroweak and the Planck scale, gauge coupling unification at high 
energy or the dark matter in the universe. However, as the supersymmetric 
partners of the SM particles have not yet been observed, supersymmetry 
must be broken at low energies, and in order to remain a viable solution 
to the hierarchy problem, this breaking must be soft, predicting massive 
superpartners around the TeV scale. Therefore, the quest for SUSY particles 
is one of the main goals of the present high-energy collider experiments such 
as the Tevatron at Fermilab or the Large Hadron Collider at CERN. 

Studying the hadron-collider phenomenology of models that go beyond 
the Standard Model requires the use of Monte Carlo event generators, based, 
on the one hand, on a proper modeling of the strong interactions to describe 
correctly the parton showering, fragmentation and hadronization, and, on 
the other hand, on the calculation of matrix-elements underlying the hard- 
scattering process. This last step requires the implementation of the complete 
set of Feynman rules associated to a given model, often one vertex at the 
time. This task, tedious if undertaken manually, has been rendered much eas- 
ier with the use of packages such as LanHep [l|, 0| or Feyn Rules [i, 0, 0|, 
starting from the Lagrangian of the theory and exporting the corresponding 
Feynman rules to one or several Monte Carlo generators. Specific tools ded- 
icated to SUSY theories, such as the Sarah package @, 0], also focusing on 
the same issue, in addition to the generation of the model mass spectrum at 
the one-loop level. However, all of these tools only address the generation of 
the Lagrangian (partially automatically or not) in the usual spacetime, and 
are not suitable for computations in superspace. 

Supersymmetric model building and the underlying calculations are often 
technically long and painstaking. The latter can be rendered easier if work- 
ing within the so-called superspace formalism jil, [9[ , a natural framework for 
SUSY model building. In this paper, we present an extension to the Feyn- 
Rules package, based on Mathematical, which provides an environment 



1 Mathematica is a registered trademark of Wolfram Research Inc. 
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to build supersymmetric theories directly in superspace. The core of our 
program is a set of basic functions allowing to expand superfields, i.e., func- 
tions defined on the superspace, in term of their component fields, the usual 
scalar, fermionic and vector fields of particle physics, in an automated way. 
As a consequence, the user has only to worry about superfield expressions, in 
general much simpler than their component fields counterparts. In addition, 
the module contains a set of predefined functions for the supercharge and 
the superderivative operators that allow the user to study the properties of 
the SUSY theory under consideration. 

Lagrangians for phenomenologically relevant supersymmetric theories can 
very often be expressed as a sum of four terms. The first two describe the 
kinetic terms for the chiral and vector supermultiplets and are independent 
of the specific model under consideration, since they are fixed entirely by 
supersymmetry and gauge invariance. For this reason, we have included into 
our code a set of functions to generate these pieces of the Lagrangian, for 
any model, in an automated way. Hence, the implementation of a SUSY the- 
ory in the FeynRules package consists only in the setting of the superfield 
content, the model parameters and the gauge symmetries of the theory as 
an input, together with the model dependent parts of the Lagrangian, i.e., 
the SUSY-breaking Lagrangian and the superpotential. This opens thus the 
way to a serious phenomenological study of entire classes of models whose 
implementation into Monte Carlo event generators was considered too com- 
plicated or too tedious so far. The FeynRules package including our su- 
perspace extension can be found, together with an up-to-date manual, at 
http:/ /feynrules. phys.ucl.ac.be 

The outline of the paper is as follows: In Section [2] we briefly review 
the main functionalities of the FeynRules package and the structure of 
the FeynRules model files. In Section [3] we introduce our implementation 
of the Grassmann and supersymmetry algebras into Mathematica, before 
turning to the more specific case of the implementation of chiral and vector 
superfields in Section HI The main part of the package, the simplification of 
superspace expressions and the generation and expansion of supersymmetric 
Lagrangians with the help of the superspace module, is discussed in Sections 
[5j M and [3 Finally, in Section [8] we illustrate the use of the package on the 
example of the implementation of the Minimal Supersymmetric Standard 
Model (MSSM) into FeynRules in terms of superfields. 
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2. The FeynRules package 



FeynRules is a Mathematica package that allows to derive Feyn- 
man rules directly from a Lagrangian |3j|. The information that the user 
needs to provide consists in the particle content and the parameters of the 
model, together with the Lagrangian that describes the interactions among 
the different particles. The Feynman rules can then be obtained auto- 
matically and the interaction vertices can be exported to various matrix- 
element generators by means of a set of translation interfaces included in 
the package. Presently, interfaces to CalcHep/CompHep loL llll 12], 
FeynArts/FormCalc (HQ], MadGraph/MadEvent fm\m O, Gj , 



Sherpa [19| and Whizard/Omega |20|, |2l| are available. In the following 
we briefly describe the basic features of the package and the model files, and 
we give a very brief introduction on how to run the code in order to derive 
the interaction vertices. For more details on both the FeynRules package 
as well as the interfaces, we refer the reader to Refs. 



The FeynRules model definition is an extension of the FeynArts 
model file format and consists in the definitions of the particles, parameters 
and gauge groups that characterize the model and the Lagrangian. Following 
the original FeynArts convention, particles are grouped into classes describ- 
ing "multiplets" having the same quantum numbers, but possibly different 
masses. Each particle class is defined in terms of a set of class properties, 
given as a Mathematica replacement list. For example, a Weyl fermion \ 
could be written as 

W[l] == { ClassName -> chi, 

Self Conjugate -> False, 

Chirality -> Left, 

Indices -> {Index [Colour] } } 

This defines a left-handed Weyl fermion (W) represented by the symbol chi. 
Note that the antiparticle is automatically declared and represented by the 
symbol chibar. The field carries an additional index labelled Colour, repre- 
senting its gauge charge under the QCD gauge group. Additional informa- 
tion, like the mass and width of the particles, as well as the U(l) quantum 
numbers carried by the fields can also be included. A complete description of 
the particle classes and properties can be found in the FeynRules manual. 

A Lagrangian is not only defined by its particle content, but also by the 
local and global symmetries defining the model. FeynRules allows to define 
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gauge group classes in a way similar to the particle classes. As an example, 
the definition of the QCD gauge group can be written 



SU3C 



== { Abelian 



-> False, 



GaugeBoson 
CouplingConstant 
StructureConstant 
Representations 



-> G, 
-> gs, 
-> f, 



-> {T, Colour} } 



where the gluon field G is defined together with the other fields during the 
particle declaration. The declaration of abelian gauge groups is analogous. 
FeynRules uses this information to construct the covariant derivative and 
field strength tensor which the user can use in the Lagrangian. 

The third main ingredient to define a model is the set of parameters which 
it depends upon. The declaration of the parameters follows the same lines 
as the declarations of the particle and gauge group classes. However, since 
the declaration of the parameters is not needed in order to understand how 
the superfield module works, we do not review it here but refer the reader to 
the FeynRules manual 

After having loaded the FeynRules package into Mathematica, the 
user can load the model and the model restrictions via the commands 

LoadModel[ filel, file2, ... ] , 

where the model can be implemented in as many files as convenient or it can 
be implemented directly in the Mathematica notebook in which case the 
list of files would be empty. The Lagrangian can now be entered directly 
into the notebook^! using standard Mathematica commands, augmented 
by some special symbols representing specific objects like Dirac matrices. As 
an example, we show the Lagrangian, 

L = -1/4 FS[G, mu, nu, a] FS[G, mu, mi, a] 
+ I chibar . sibar[mu] . DC[chi, mu] , 

where FS [G , mu, nu , a] and DC [chi , mu] denote the SU(3)c field strength 
tensors and covariant derivatives automatically defined by FeynRules, re- 
spectively. At this stage, the user can perform a set of basic checks on the 



2 Alternatively, the Lagrangian can also be included in the model file, in which case it 
is directly loaded together with the model file. 
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Lagrangian (hermiticity, normalization of kinetic terms, . . . ), or directly pro- 
ceed to the derivation of the Feynman rules via the command 

verts = FeynmanRules [ L ] . 

FeynRules then computes all the interaction vertices associated with the 
Lagrangian L and stores them in the variable verts. The vertices can be used 
for further computations within Mathematica, or they can be exported to 
one of the various matrix-element generators for further phenomenological 
studies of the model. The translation interfaces can be directly called from 
within the notebook, e.g., for the FeynArts interface, 

WriteFeynArtsOutput [ L ] , 

This will produce a file formatted for use in FeynArts. All other interfaces 
are called in a similar way. 

The FeynmanRules function, as well as the interfaces that call it, require 
the Lagrangian to be written in four- dimensional space time. For supersym- 
metric theories, however, the most natural and most convenient way to write 
a Lagrangian is in terms of superfields, i.e., in terms of the explicit super- 
multiplets of the supersymmetry algebra. For this reason, we have extended 
the FeynRules package by a superspace module, which allows to write a 
Lagrangian in terms of superfields and to transform the latter into a four- 
dimensional Lagrangian involving component fields only, hence expressing 
the superfield Lagrangian in a way that can be processed immediately by the 
FeynmanRules function. This new module will be described in detail in the 
next sections. 



3. Superspace and supersymmetry algebra 

3.1. Grassmann variables and fermionic fields 

Supersymmetric theories are naturally formulated in superspace [§, [9|, 
an extension of the ordinary spacetime, defined, for N = 1 supersymmetric 
theories by adjoining a Majorana spinor (9 a , 9 a ) to the usual spacetime coor- 
dinates x^, where 9 and 9 are (Grassmannian) two-component Weyl fermions. 
Using the recent implementation of Weyl fermions into FeynRules 22j], the 



9 variables are hardcoded in the superspace module following the example of 
Section [2], as 
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W[xlOOO] == { 

TeX -> \[Theta] , 

ClassName -> theta, 

Chirality -> Left, 

Self Conjugate -> False, 

Unphysical -> True} 

This field and its right-handed counterpart 9 can be accessed, after loading 
FEYNRULES in MATHEMATICA, through the symbols theta and thetabar. 
By convention, the spin indices are assumed to be lowered for both the 9 and 
9 variables, 

theta [alpha] 6 a and thetabar [alphadot] 6& , 

where we relate undotted and dotted indices to left and right-handed fermions. 
Spin indices can be raised and lowered using the rank-two antisymmetric ten- 
sors e aj 3 and 

ff* = e a %, 9 a =e aP 9^ 9 & = e^fy and 9 & = e-^, (1) 

with = 1, e 12 = —1, — 1 an d e 12 = — 1. This strictly follows the 
conventions of Ref. 24J, and so does the complete superspace module of 
FeynRules. Regardless of the left or right-handed nature of the spin indices, 
the tensors with lower and upper indices are implemented into FeynRules 
as Deps and Ueps, respectively. An example is in order, 

Ueps [alphadot, betadot] thetabar [betadot] <H> £ Q/3 #g = 9 a , 
Deps [alpha, beta] Ueps [beta, gam] theta[gam] -H- e a pe l3l 9 1 = 9 a . 

The Levi-Civita tensors allow us to fix our conventions with respect to the 
summation on the spin indices, and we define the dot product of two Weyl 
spinors as 

A • A' = X a X' a = e^Xp\' a , and x • x' = XaX' & = e^XaX'p , (2) 

where we have introduced some generic left and right-handed spinors A^ and 
X^', respectively. 

Computations in superspace require to keep track not only of the position 
of the different spin indices, but also of the ordering of the fermions. The 
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FeynRules superspace module addresses this issue in the following way. 
First, any explicit spin index will be considered, by default, as lowered. Sec- 
ond, we have implemented the environment nc [chain] , where chain stands 
for an ordered sequence of fermions with lower spin indices. As simple ex- 
amples, scalar products such as in Eq. (jSJ) are implemented as 

-H- nc [lambda [a] , lambdaprime [b] ] Ueps[b,a] , 

a-_ ( 3 ) 
X ■ X = £ XbXol nc [chi [bd] , chiprime [ad] ] Ueps[bd,ad] . 

3.2. Supercharges and superderivatives 

In N = 1 supersymmetric theories, there is a single spinorial SUSY gen- 
erator, transforming as a Majorana spinor {Q a ,Q a ) under Lorentz transfor- 
mations. Any translation in superspace can then be parameterized as 

G(x,6,d) = e i \- x ^ +e - Q+Q ^ , (4) 

where P M denotes the generator of space-time translations. A generic super- 
symmetric transformation then corresponds to G(0,e,e), where (e a ,e a ) de- 
notes a (spinorial) transformation parameter. Applying this transformation 
to a generic superspace point, G(0,e,e)G(x,9,9), we can calculate the vari- 
ations of the spacetime coordinates x and of the Grassmann variables 9 and 
9 using the Baker-Campbell-Hausdorff formula and the (anti) commutation 
relations among the generators. Comparing with a direct application of the 
supersymmetric generator z(e • Q + Q ■ e) on the coordinates x, 9 and 9, 
we obtain the action of the generators as differential operators acting on 
functions on superspace. Similarly, starting from an action from the right, 
G(x, 9, 9)G(0, e, e), we obtain the expressions of the superderivatives D a and 



D a . In our conventions 24j, these four quantities are given by 



(5) 



The matrices a 11 , and for later reference a 11 , are the four- vectors correspond- 
ing to the usual Pauli matrices. Note that, since the actions from the left 
and from the right commute, the supercharges and the superderivatives must 
anticommute. 
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The generators Q a and Qa of the supersymmetric transformations, as 
well as the superderivatives D a and D a , are called in FeynRules via the 
functions, 

QSUSY [expression, alpha] <H> Q a (expression) , 
QSUSYBar [expression, alphadot] -H- Q a (expression) , 
DSUSY [expression, alpha] -H- D a (expression) , 
DSUSYBar [expression, alphadot] -H- D a (expression) . 

It is assumed that all the indices appearing in expression must be written 
explicitly, using the environment nc presented above in order to keep track of 
the correct fermion ordering. The only exception corresponds to the case of 
a single fermion, the ordering of the fermions being trivially irrelevant, and 
the nc environment can be omitted. Hence, both 

DSUSY [chi [alpha] , beta] and DSUSY [nc [chi [alpha] ] , beta] 

give a correct answer, whilst for two fermions, the correct synthax is 
DSUSY [nc [chi [alpha] , lambda [beta] , gamma] , 

and not 

DSUSY [chi [alpha] lambda [beta] , gamma] . 

4. Implementing superfields into FeynRules 

4-1. Generic superfields 

Any function <&(x, 9, 9) defined on N = 1 superspace is called a superfield 
and can be expanded into a Taylor series with respect to the anticommuting 
coordinates 9 and 9. Since the square of an anticommuting object vanishes, 
the series has only a finite number of terms and the most general expression 
for a scalar superfield i^l 

9, 9) = z(x) + e-£(x) + 9-C(x) + 9-9 f{x) + 9-9g{x) 

+ 9a»9V fl (x) + 9-9 9-u{x) + 9-9 9-p(x) + 6-6 9-9d(x) , ^ 



3 The extension to non-scalar superfields is immediate. 
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Table [D Superspace functionalities 

theta[al] The Grassmann variable 9 a , where al denotes the 

carried spin index. 

thetabar [aldot] The Grassmann variable 9 a , where aldot denotes 
the carried spin index. 

Ueps[i,j] The rank- two antisymmetric tensor e 1 - 7 . 

Deps[i,j] The rank- two antisymmetric tensor e^-. 

nc[seq] Ordered sequence of fermionic field(s), labelled by 

seq, where each field carries explicitly its indices. 

QSUSY[exp,al] Calculates the action of the supercharge Q a on 
the expression exp, where all indices must be writ- 
ten explicitly, using the nc environment if relevant. 
The symbol al is related to the spin index of the 
supercharge. 

QSUSYBar [exp, ad] Calculates the action of the supercharge Q a on 
the expression exp, where all indices must be writ- 
ten explicitly, using the nc environment if relevant. 
The symbol ad is related to the spin index of the 
supercharge. 

DSUSY[exp,al] Same as QSUSY, but for the super derivative D a . 
DSUSYBar [exp , ad] Same as QSUSYBar, but for the superderivative D a . 
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where the various coefficients of the expansion form a so-called supermultiplet 
and are referred to as the component fields of the superfield. They correspond 
to the usual scalar, fermionic and vector degrees of freedom used in particle 
physics, z, f, g and d denote complex scalar fields, whilst £, £, uj and p 
denote complex Weyl fermions and is a complex vector field, leaving a 
total of 16 bosonic and 16 fermionic degrees of freedom. 

Using the nc environment, the generic scalar superfield of Eq. (J5]) can be 
implemented into FeynRules in a straightforward way as 

Phi = 

z + <H> z 

nc [theta[sp] , xi[sp2]] Ueps[sp2,sp] + -H- 9-£ 

nc [thetabar [spd] , zetabar [spd2] ] 

Ueps[spd,spd2] + <H- 9-Q 

nc [theta [sp] , theta[sp2]] Ueps[sp2,sp] f + <H- 9 -Of 
nc [thetabar [spd] , thetabar [spd2] ] 

Ueps [spd, spd2] g + ^9-9g 
nc [theta [sp] , thetabar [spd] ] Ueps[sp2,sp] 

Ueps [spd2, spd] si [mu , sp2 , spd2] V[mu] + <-> 9a fi 9V fl (?) 
nc [thetabar [spd] , thetabar [spd2] ] 

Ueps [spd, spd2] Ueps[sp2,sp] 

nc [theta [sp] , omega [sp2]] + ^-9-9 9-uo 

nc [theta [sp] , theta [sp2]] Ueps[sp2,sp] 
nc [thetabar [spd] , rhobar [spd2] ] 

Ueps [spd, spd2] + ^9-9 9- p 

nc [theta [sp] , theta [sp2]] Ueps[sp2,sp] 
nc [thetabar [spd] , thetabar [spd2] ] 

Ueps [spd, spd2] d ^9-99-9d, 

where we assume that the component fields have been defined properly in 
the FeynRules model file (See Section [2]). In Eq. (J7|) we have made use 
of the object si [mu, alpha, alphadot] , corresponding to the Pauli matrices 
cr^aa- Note that the appearance of the Levi-Civita tensor Ueps is due to our 
convention that all spin indices carried by fermionic fields are lowered. Sim- 
ilarly, we could have used the matrices sibar [mu, alphadot , alpha] defined 
with two upper spin indices, a ma since 

e^a» a ^9 = -h<?^ a 9 a . (8) 
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(9) 



The second option would have allowed to get an expression free from any 
e-tensor. 

Even though Eq. (J7J) is the canonical form in which superfields are rep- 
resented inside the code, the form of Eq. ([7]) can however be very painful 
to use in practice. For this reason, the package contains a function denoted 
ncc, which has the same effect as nc, but all spin indices and e-tensors are 
suppressed. The code then assumes that, within a given ncc environment, 
all the suppressed indices are contracted according to the convention of Eq. 
03]), and outputs the result into the canonical form. As an example, using 
the ncc environment, Eq. (J7|) can be written in the more compact form, 

Phi = 

z + <-»■ z 

ncc[theta, xi] + -H> 

ncc [thetabar , zetabar] + <H> 

ncc[theta, theta] f + <H- 9-9 f 

ncc [thetabar , thetabar] g + <H-9-9g 

ncc [theta, si[mu], thetabar] V[mu] + -H- Ou^OV^ 

ncc [thetabar , thetabar] ncc [theta, omega] + ^9-9 9-u) 

ncc[theta, theta] ncc [thetabar , rhobar] + ^-9-9 9-p 

ncc[theta, theta] ncc [thetabar , thetabar] d ^-9-9 9-9d. 

The number of degrees of freedom of a generic superfield is too large to 
match those of the supermultiplets representing the N = 1 supersymmetric 
algebra. Two special cases, with fewer degrees of freedom, are in general 
enough to build most of the phenomenologically relevant supersymmetric 
theories. These so-called chiral and vector superfields will be discussed in 
the next sections. 

4-2. Chiral superfields 

Left and right-handed chiral superfields are superfields $^ and that 
satisfy the constraints 

D & <$> L (x, 9,9) = and D a ® R (x, 9, 9) = , (10) 

where the superderivatives have been introduced in Eq. (J5J). Since the su- 
perderivatives anticommute with the supercharges, the constraints are pre- 
served by supersymmetry transformations. The most general solutions to 
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the constraints can be written as 

9) = <f>{y) + y/29 ■ j>{y) - 9 ■ 9F(y) , 
§ R {y\ 9) = 0(yt) + . ^( y t) _ . F ( y t) , 

where we have introduced the varia y* 1 = x^ — i6a^9 and its hermitian 
conjugate (y^Y = x^ + iOa^d. Chiral superfields are appropriate to describe 
the N = 1 matter supermultiplets containing one fermion and one scalar 
field. Indeed, the physical degrees of freedom of $l and consist of a 
single complex scalar field and a single two-component fermion ip, the chi- 
rality of the Weyl fermion providing the name for this kind of superfield. 
The additional scalar field F has a mass dimension [F] = 2 and does hence 
not correspond to a physical degree of freedom, but it is necessary to restore 
the equality between the numbers of fermionic and bosonic degrees of free- 
dom off-shell. Furthermore, these so-called F-terms can be used to break 
supersymmetry spontaneously if they develop a vacuum expectation value. 
As we will discuss in Section 17.2^ these fields are non-propagating and can 
be eliminated through their equations of motion. 

Left and right-handed chiral superfields can be declared in FeynRules 
following in the FeynRules model file in a way similar to the ordinary fields, 
e.g., 

M$Superf ields 
CSF[1] == { 
ClassName 
Chirality 
Weyl 
Scalar 
Auxiliary 

CSF[2] == { 
ClassName 
Chirality 
Weyl 
Scalar 

} • 



4 The \pl factors and the minus signs are purely conventional. 
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= { 

-> PHI, 
-> Left, 
-> psi, 
-> z, 
-> FF}, 



-> OMEGA, 
-> Right, 
-> xibar, 
-> zz} 



These Mathematica commands declare two chiral superfields, the label of 
the particle class CSF referring to chiral superfields. We consider the example 
of both a left and right-handed chiral superfield called PHI and OMEGA. The 
fermionic degrees of freedom of the superfields are psi and xibar, while the 
scalar degrees of freedom are denoted by z and zz. In addition to all the 
options available to all particle classes and briefly reviewed in Section 121 such 
as QuantumNumbers and Indices, the user must set the Chirality option 
to Left or Right, related to the chirality of the superfield. Moreover, it is 
mandatory to link a chiral superfield to its fermionic and scalar components 
through the options Weyl and Scalar, respectively. Note that the component 
fields have to be declared independently in the FeynRules model file (see 
Section[2]for a brief review or Refs. ji, 22] for a detailed description). We only 



emphasize here that the options for the component field classes, like Indices , 
QuantumNumber, etc., must be identical to the corresponding options of the 
superfield. The only exception to this rule are the F-terms, where the user 
can either point to an already declared complex scalar field via the Auxiliary 
options (as for CSF[1] in the example above, where the auxiliary field is 
denoted FF), or he leaves this option unspecified and FeynRules creates 
internally a symbol for the auxiliary field (as for CSF [2] in the example 
above). All the allowed options for the declaration of a chiral superfield are 
summarized in Table |2j 

4-3. Vector superfields in the Wess-Zumino gauge 

Besides matter fields, gauge theories contain also real vector fields describ- 
ing the gauge bosons. Since the chiral superfields defined in the previous sec- 
tion do not contain any vector degree of freedom, they cannot be sufficient to 
describe supersymmetric gauge theories. We therefore introduce in this sec- 
tion the so-called vector (or gauge) supermultiplets, i.e., the representations 
of the N = 1 SUSY algebra containing one massless gauge boson together 
with the corresponding fermionic degree of freedom. These multiplets can be 
described by vector superfields, defined by the reality condition, 

V = V* . (12) 

The constraint (Tl2|) on its own is not enough to reduce the number of de- 
grees of freedom of the generic superfield of Eq. to the required number 
and leads to a proliferation of unphysical fields that can be eliminated by a 
suitable gauge choice. A convenient choice is the so-called Wess and Zumino 
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Table [21 Chiral superfield class options 


ClassName 


Defines the symbol by which a class is represented. 


Chirality 


Defines the chirality, Left or Right, of the chiral 
superfield. 


Weyl 


Contains the symbol of the fermionic component of 
the chiral superfield. 


Scalar 


Contains the symbol of the scalar component of the 
chiral superfield. 


Optional attributes 


Auxiliary 


Contains the symbol of the auxiliary component of 
the chiral superfield. If absent, FeynRules gen- 
erates the F-term automatically. 


Indices 


The list of indices, different from Lorentz and spin 
indices, carried by the superfield and all its com- 
ponent fields. 


QuantumNumbers A replacement rule list, containing the U(l) quan- 
tum numbers carried by the class. 
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gauge, in which a vector superfield is expressed as 

y wz = Oo^Ov^ + i9 - 9 9 -X-i9 - 9 9 ■ X+^9 - 9 9 -9D . (13) 

The component fields of a vector superfield are hence a real vector field 
and a Majorana fermion (A a , X a ). Similar to the case of the chiral superfield, 
the multiplet also contains a non-propagating auxiliary scalar field D with 
mass dimension [D] = 2 that is necessary to restore the equality between 
the numbers of fermionic and bosonic degrees of freedom off-shell. These so- 
called .D-terms can again be eliminated through their equations of motion. 

Vector superfields in the Wess and Zumino gauge can be implemented 
into FeynRules in the same way as chiral superfield. Let us illustrate this 
by an example, 

M$Superf ields = { 
VSF[1] == { 

ClassName -> PHIV, 

GaugeBoson -> X, 

Gaugino -> lambda} 

} . 

In this example, PHIV denotes a vector superfield, whose vector component 
is denoted by X and the fermionic component by lambda. The label of the 
particle class VSF refers to its vector superfield nature. In addition to the 
usual options for particle classes, the user must associate to a vector superfield 
its bosonic and fermionic components by setting the options GaugeBoson and 
Gaugino. Note that the corresponding vector field and the Weyl fermion 
must be declared independently in the FeynRules model file. As in the 
case of chiral superfields, the D-terms can either be defined explicitly via 
the Auxiliary option, or this option can be omitted and FeynRules will 
declare the D-terms internally. All the allowed options for the declaration 
of a vector superfield are summarized in Table [31 In addition, each vector 
superfield can be linked to a gauge group through the option Superfield 
which has been added to the gauge group class in a similar way as the we 
assigned the gluon field G to the QCD gauge group SU3C in Section [2J As an 
example, we could associate the superfield PHIV above to an abelian gauge 
group called U1X by declaring the gauge group as 

U1X == { 
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Table [3t Vector superfield class options 


ClassName 


Defines the symbol by which a class is represented. 


GaugeBoson 


Contains the symbol of the vector field associated 
to the vector superfield. 


Gaugino 


Contains the symbol of the gaugino component of 
the vector superfield. 


Optional attributes 


Auxiliary 


Contains the symbol of the auxiliary component of 
the superfield. If absent, FeynRules generates 
the D-term automatically. 


Indices 


The list of indices, different from Lorentz and spin 
indices, carried by the superfield and all its com- 
ponent fields. 


New gauge g 


roup options 


Superfield 


This option points to the ClassName of the vec- 
tor superfield associated to the gauge group. If 
the Superfield option is present, the option 
GaugeBoson becomes optional. If both options are 
there, they must be consistent. 


Abelian -> True, 
CouplingConstant -> gX, 
Superfield -> PHIV} . 



5. Simplification of superspace expressions 

5.1. Simplification of expressions involving Grassmann numbers 

In the previous sections we have discussed the implementation of super- 
fields into FeynRules. The user can either use the predefined classes CSF 
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and VSF for chiral and vector superfields, or implement superfields directly in 
terms of the component fields, as in the example of Eq. (0). In that case, how- 
ever, the required synthax for the implementation can be complex, and as a 
consequence the Mathematica output could be difficult to read. Similarly, 
many calculations such as those involving the supersymmetric generators or 
the superderivatives could lead to rather long expressions after introducing 
the second rank antisymmetric tensors and the nc environment, which are 
necessary in order to have the fermions ordered correctly and all the spin in- 
dices carried by fields lowered, as is the convention in the superspace module 
of FeynRules. To bypass this issue, it might be useful to form dot products 
of spinors, and more specifically those involving the Grassmann variables 9 
and 9. This can be achieved with the help of the ToGrassmannBasis function, 
which proceeds in two steps. First, dot products are formed and simplified 
using relations among the Grassmann variables, such as 

e a 9 p = —9 ■ 9e afi , 9 & 9 $ = -9 ■ 9e^ , 9 a 9 & = -9a»9a^ a , ... (14) 
2 2 2 M K J 

After this step, any function of the component fields is written in terms of 
a restricted set of scalar products involving Grassmann variables and Pauli 
matrices, forming hence a basis in which any superfield expression can be 
expanded. 

Even though this procedure in principle solves the problem of simplifying 
superfield expressions, we have to deal at this stage with a purely technical 
issue. Expressing everything in terms of the basic objects results in Math- 
ematica expressions that are equal up to the names of contracted indices, 
e.g., 

Dot [theta[al] , theta[al]] - Dot [theta[be] , theta[be]] . (15) 

Although this difference is manifestly zero, the two terms represent different 
patterns in Mathematica, and hence the cancellation does not take place. 
In general, such a situation can occur relatively often with the ToGrassmann- 
Basis function, and therefore the internal index naming scheme is optimized 
by the ToGrassmannBasis function after all the dot products are formed. 
Similar terms are collected and summed, and hence the readability of the 
final results is highly improved. The simplification function presented above 
can be called in FeynRules through the command 

ToGrassmannBasis [expression] , 
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where expression is any function of the component fields. The ToGrass- 
mannBasis function then expresses expression in the basis defined by the 
objects in Eq. (TH|) . 

As an example, the function ToGrassmannBasis could significantly help 
to improve the readability of the most general series in the 9 and 6 variables 
presented in Eq. (J7J). Using the Mathematica variable Phi defined in Eq. 
dZj), the command 

ToGrassmannBasis [Phi] 

allows the user to obtain an expression much closer to the original form of 
Eq. ©, 

z + theta[sp] .xi [sp] + zetabar [spd] .thetabar [spd] + 
f*theta[sp] . theta[sp]+g*thetabar [spd] .thetabar [spd] + 
theta[sp] .thetabar [spd] *si [mu, sp, spd] *V[mu] + 
rhobar[spd] .thetabar [spd] *theta[sp] .theta[sp] + 
theta[sp] . omega [sp] *thetabar [spd] .thetabar [spd] + 
d*theta[sp] . theta[sp] *thetabar [spd] . thetabar [spd] . 

It is important to note that the ToGrassmannBasis function can also work 
on spinorial or tensorial expressions, i.e., expressions with uncontracted spin 
(or vectorial) indices. After an application of the simplification module, the 
(upper or lower) free index could be attached to a single fermion, or to a 
chain containing one fermion and a given number of Pauli matrices, 

X a , X^aa ^ (XO • , (T^^XP ^ {^X) a , ■■■ , (17) 

where, in the examples above, x denotes a generic left-handed Weyl fermion. 
The ToGrassmannBasis function has been implemented so that those chains 
are formed and stored in the TensDot2 environment following the pattern 

TensDot2[ chain ] [pos, chir, name] , 

where chain is a sequence of one Weyl fermion and possibly one or several 
Pauli matrices, pos is the up or down position of the free spin index, chir its 
chirality, i.e., its dotted or undotted nature, and name its name. The three 
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examples of Eq. (1171) could be implemented as 

X a -H- ToGrassmannBasis [ nc[chi[b]] * Ueps[a,b] ] , 

{x a>1 )a ToGrassmannBasis [ nc[chi[b]] * si[mu,a,ad] * 
Ueps[a,b] ] , 

[a^a v x) ToGrassmannBasis [ nc[chi[b]] * s i [mu , a , ad] * 
sibar [nu, ad,b] ] . 

We recall that the nc environment is mandatory as soon as we are dealing 
with fermions and that by convention all spin indices carried by fermionic 
fields are considered to be lowered. We then obtain 

nc[ TensDot2 [chi [a] ] [up,Lef t , a] ] , 
nc [ TensDot2 [chi [a] , si [mu, a, ad] ] [down, Right , ad] ] , 
nc [ TensDot2 [si [mu,a,ad] , sibar [nu, ad, b] , 
chi [b]] [down, Left, a] ] . 

To conclude this section, let us comment on the index optimization rou- 
tine used by the ToGrassmannBasis function, because there might be times 
where the user wants to handle the optimization of the index naming scheme 
without forming scalar products involving Grassmann variables, i.e., without 
calling the ToGrassmannBasis function. The standalone version of the index 
optimization that consistently renames the indices of an expression is called 
as 

Optimizelndex [expression, list] , 

where list is an optional list of variables carrying indices to be included into 
the index renaming and that are neither a field nor a parameter included 
in the global variable of FeynRules M$Parameters. As an example, the 
application of the optimization of the index naming scheme on the sum of 
scalar products of Eq. ( II 5p can be performed with the command 

Optimizelndex [ Dot [theta[al] ,theta[al] ] - 
Dot [theta [be] , theta [be] ] ] , 

which simplifies to zero. Equal terms have now been consistently subtracted 
and are not repeated anymore. 



21 



5.2. Simplified input format for superfield expressions 

In Section [3] we introduced a canonical form in which every superspace 
function can be expressed. This canonical form is defined by the two simple 
rules 

1. all spin indices carried by fermion fields are considered lowered, 

2. the ordering of the fermions is implemented via the nc environment. 

As we already discussed, this canonical form usually requires the explicit 
use of the e-tensors, and for this reason, we have introduced in the pre- 
vious section the ToGrassmannBasis function which reduces any superfield 
expression to the basis of Grassmann variable monomials defined by Eq. (|14p . 
Conversely, this basis can also be used to input superspace expressions in an 
easier way. The rules for this input format are discussed in the rest of this 
section. 

First, the dot products of spinors, connected or not with the help of Pauli 
matrices, are always written as 

f erml [spl] . f erm2 [sp2] chain [spl , sp2] , 

where the symbols ferml and ferm2 denote the two fermions and chain 
contains a series of Pauli matrices linking the two spin indices spl and sp2. 
As a simple example, a possible connected product of two left-handed spinors 
X and ip would be 

XC^a^tp -H- chi [spl] .psi [sp2] si [mu, spl , spd] sibar [mi, spd, sp2] , 

the symbol chain being here equal to si [mu, spl , spd] sibar [nu, spd, sp2] . 

In the cases where there is no chain of Pauli matrices present in the expres- 
sion, the spin indices carried by the fermions must be equal, the dot product 
being hence a regular scalar product of fermions, as those introduced in Eq. 
(T5]). A few additional examples can be found among the different terms of 
Eq. (|T6|) . Second, any fermionic expression carrying a free spin index must 
use both the nc environment as well as the TensDot2 structure, following the 
synthax explained in Section 15.11 

Even though this format allows to input superspace expressions in an 
easier way, most of the functionalities of the superspace module require the 
input expressions to be given in the canonical form of Section |3j In the next 
section we therefore introduce a function that allows to convert an expression 
written in terms of the basis objects of Eq. ( TPfj) to its canonical form. 
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5.3. Reverting simplifications: back to the nc environment 

One drawback of the optimization routines is that several of the func- 
tions included in the superspace module of FeynRules, such as the QSUSY 
or DSUSY routines presented in Section [31 require expressions including the 
nc environment, second rank antisymmetric tensors, lower spin indices for 
fermions and not any explicit scalar products. A simple solution is provided 
by the Tone function. Indeed, if one needs to perform additional operations 
on simplified expression within the superspace, it is recommended to first use 
the Tone function, which allows to transform simplified expressions back to 
their original form in terms of the nc environment and the epsilon tensors. 
Hence, Tone [Dot [theta [a] , theta [a] ] ] would lead to 



where the second spin index has been automatically generated by Feyn- 
Rules. Sometimes, it might be useful to call the Optimizelndex function 
right after the use of the Tone module in order to get shorter and simpler 
expressions. 

Note that the Tone function allows us at the same time to enter expres- 
sions directly with the use of the simplified synthax presented in Sections 15. II 
and 15.21 and to convert them to their canonical form. 

6. Manipulating superfield expressions 

6.1. From superfields to Lagrangians 

The main advantage of writing down supersymmetric Lagrangians in 
terms of (chiral and vector) superfields rather than in terms of the com- 
ponent fields is the size of the corresponding expressions. This is illustrated 
in the following simple example. Let us consider a left-handed chiral super- 
field $ whose scalar, fermionic and auxiliary components are denoted by </>, 
ip and F . Performing the scries expansion of $t$ in terms of the Grassmann 
variables, it can be shown that the coefficient with the highest power in 9 
and 9, i.e., the coefficient of the 9 ■ 99 -9 term, is SUSY invariant and hence 
a good candidate for a Lagrangian density describing the free component 
fields. Indeed, we get 



nc [theta [a] , theta [sp$l] ] Ueps [sp$l , a] 



£ = $ f <I> 




(18) 
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Table [4} Simplification functions 

ToGrassmannBasis [exp] 

This function allows to express any function exp 
of the component fields in the basis of Grassmann 
monomial defined by Eq. f fT4|) . 

Optimize Index [exp , list] 

This function allows for the optimization of the in- 
dex naming scheme used in the expression exp. The 
optional argument list is a list of variables carry- 
ing indices, which are neither a field nor included 
in M$Parameters. 

Tone [exp] This function transforms simplified expressions 

back to their original form, with lower spin indices, 
epsilon tensors, etc... 

New environments 

TensDot2 [chain] [pos,chir,name] 

This environment contains a sequence, labelled by 
chain, of one Weyl fermion and possibly one or 
several Pauli matrices. The symbols pos, chir and 
name are the up or down position, the chirality and 
the name of the free index. 
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The expansion of a superfield polynomial into a series in the Grass- 
mann variables can be automatically performed in FeynRules via the SF2- 
Components function. Since this function allows in the same way to transform 
superfield expressions into four-dimensional Lagrangians, it is one of the most 
important functions available and at the heart of the superspace module of 
FeynRules. For any given function of chiral and vector superfields denoted 
expression, the correct syntax to use is simply 

SF2Components [ expression ] . 

The SF2Components function expands all the superfields appearing in ex- 
pression in terms of their component fields and the x M spacetime coordinates 
(rather than the variable related to chiral superfields). In a second step, 
scalar products of Grassmann variables are simplified and the expression 
is reduced to the basis defined by Eq. (I14p using the ToGrassmannBasis 
function. During this procedure representation matrices of the Lie algebra of 
the gauge groups are introduced and simplified, using, e.g., the commutation 
relations between the generators. The output of the SF2Components function 
consists in a list of two elements, 

{ Full series , List of the nine coefficients } . (19) 

The first element of this list, labelled here Full series, is the full series 
expansion in the Grassmann variables which could also have been directly 
obtained with the GrassmannExpand function, 

GrassmannExpand[ expression ] 

The second element of the list of Eq. ( lT9"j) is a list containing the nine coef- 
ficients of the series. The first part of this last list is the scalar piece of the 
series, i.e. the terms independent of the 9 and 9 variables, while the other 
elements are the coefficients of the 9 a , 9a fJ- 9, 9 ■ 9, 9 ■ 9, 9 ■ 99^, 9 ■ 99 a and 
9-99-9 terms, following this ordering. Each of them can also be obtained 
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using the functions, 



ScalarComponent [ expression ] , 
ThetaComponent [ expression ] , 
ThetabarComponent [ expression ] , 
ThetaThetabarComponent [ expression ] , 
Theta2Component [ expression ] , (20) 
Thetabar2Component [ expression ] , 
Theta2ThetabarComponent [ expression ] , 
Thetabar2ThetaComponent [ expression ] , 
Theta2Thetabar2Component [ expression ] . 

The results obtained by each of these functions are thus all independent of 
the Grassmann variables. It is important to note that each of these functions 
calls the SF2Components module, the complete series being thus recalculated 
at each function call. Hence, if several of the coefficients have to be computed, 
it is much faster to store and re- use the results of the SF2Components function 
than to call all the XXXXComponent functions individually. Furthermore, these 
functions can also be used on expressions in terms of the component fields. 
In this case, we recall that the use of the Tone environment is mandatory. All 
the available functions to reexpress superfield expressions in terms of their 
component fields are summarized in Table 

As an example, let us derive the Lagrangian of Eq. (fT5|) with FEYN- 
Rules, starting from the superfield expression, recalling that the superfield 
PHI has been defined in Section I4~2l This can be done via the Mathematica 
command 

Theta2Thetabar2Component [PHIbar PHI] , 

and one obtains indeed the correct free Lagrangian for the component fields, 

(del[z, mu] *del [zbar , mu])/2 - 
(zbar*del [del [z, mu] , mu])/4 - 

(z*del [del [zbar, mu] , mu])/4 + FF*FFbar - (21) 
(1/2) *del [xi [sp] ,mu] . xibar [spd] * si [mu,sp,spd] + 
(1/2) *xi [sp] .del [xibar [spd] ,mu] * si [mu,sp,spd] , 

where del stands for the spacetime derivative. 
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Table [5 


\ From superfields to particles 


SF2Components [exp] 


Expands the superfield expression exp in terms 




of its component fields and simplifies the prod- 




ucts of Grassmann variables. The result is a list 




of two elements, the complete series and a new 




list with all the individual coefficients. 


Shortcuts to the individual component fields 


The full series. 


GrassmannExpand [exp] 


The scalar term. 


ScalarComponent [exp] 


The 9 term. 


ThetaComponent [exp] 


The 9 term. 


ThetabarComponent [exp] 


The 9u9 term. 


ThetaThetabarComponent [exp] 


The 9 2 term. 


Theta2Component [exp] 


The 9 2 term. 


Thetabar2Component [exp] 


The 9 2 9 term. 


Theta2ThetabarComponent [exp] 


The 9 2 9 term. 


Thetabar2ThetaComponent [exp] 


The 9 2 9 2 term. 


Theta2Thetabar2Component [exp] 
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6.2. Super symmetric transformation laws 

The supersymmetric transformation laws of a superfield can be obtained 
by using the explicit representation of the supercharges given in Eq. (J5J). 
Hence, considering an infinitesimal supersymmetric transformation of a (spino- 
rial) parameter (e a , e°), a superfield transforms as 

$^$ + <5 e $ = $ + z(e-Q + Q-e)$ . (22) 

After expanding the superfield in terms of its component fields and using 
Eq. ()5]), we can immediately read off the transformations of the component 
fields. The operator 5 t is implemented into the superfield module via the 
DeltaSUSY function, 

DeltaSUSY [ expression , epsilon ] , 

where expression can be any function of superfields and/or component fields 
and epsilon refers to the supersymmetric transformation parameter, with- 
out any spin index. There are ten such parameters predefined in the su- 
perfield module, labelled by epsx with x being an integer between zero and 
nine, which the user can use at his convenience. This parameter being a Ma- 
jorana fermion, it is enough to only provide the associated two-component 
left-handed spinoJE The output of the DeltaSUSY module is the full series 
expansion in the Grassmann variables. Finally, the functions given in Eq. 
(120]) allow for the identification of the variations of the various component 
fields. 

Let us illustrate the use of this function on the explicit example of chiral 
superfields, and let us start with the superfield PHI introduced in Section 
14.21 We then calculate its variation under a supersymmetric transformation 
of parameter ei, 

DeltaPHI = DeltaSUSY [ PHI , epsl ] , 

where we have used the supersymmetric transformation parameter epsl. 
The scalar, 9 a and 9 ■ 9 coefficients of DeltaPHI are finally identified with the 
variations of the scalar component 5 ei z, the fermionic component 5 tl ip and of 



5 The objects epsx are in fact implemented into FeynRules as left-handed Weyl 
fermions. The corresponding right-handed objects are also available under the name 
epsxbar. 
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the auxiliary component 8 ei F of the superfield. The individual components 
can be easily extracted using the shortcuts of Eq. ( I2(jp . and one recovers the 
well-known textbook expressions, 



8 ei z 
5 ei ijj 
5 ei F 



V 7 ^! ■ if, , 



(23) 



where z, ip and F are the scalar, fermionic and auxiliary component of the 
considered superfield. 

7. Implementing supersymmetric Lagrangians into FeynRules 

7.1. Supersymmetric Lagrangians 

In this section we describe how to generate in an automated way (parts 
of) the Lagrangian describing the interactions between the different compo- 
nent fields of chiral and vector supermultiplets. The first part of this section 
is devoted to a brief review on how to construct supersymmetric Lagrangians. 
The most general supersymmetry-conserving Lagrangian describing the in- 
teractions between the various multiplets can be written as a sum of three 
pieces, 



where £ c hi ra i an d ^Yang-Miiis contain the kinetic terms as well as the gauge 
interactions of the different particles, and £ SBper w is the superpotential de- 
scribing the interactions between the different chiral supermultiplets. Note 
that, since supersymmetric particles have not yet been observed, supersym- 
metry must be broken at low energies, which renders the superpartners heavy 
in comparison to their Standard Model counterparts. The corresponding La- 
grangian does not involve superfields, but only some of the component fields, 
and so we exclude it from the present discussion and only concentrate on 
Lagrangians where supersymmetry is unbroken. 

The first Lagrangian in Eq. (124)) . £ c hirai> contains the kinetic terms as well 
as the gauge interactions of the chiral superfields. It is completely fixed by 
gauge invariance and reads 




(24) 



chiral 



(25) 



e 



66 



6 
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where a sum over all chiral superfields $ l and vector superfields = V^ a T a 
is understood. Furthermore, T a denote the representation matrices of the 
supermultiplet and gj is the associated gauge coupling constant. It can 
be shown that the 9 -99 -9 coefficient of the expansion in the Grassmann 
variables of the expression inside the squared brackets is invariant under 
supersymmetry transformations and reproduces for gauge singlet superfields 
the Lagrangian (|18[) for a free chiral supermultiplet. This expression is hence 
a good candidate for a Lagrangian density. In general, if the fields carry 
some gauge charges, the expansion of Eq. ( |25|) yields 



C 



chiral 



D^\D^ + F-F l 



iV2gf\ ja ■ ^T a l - iV2g j( f>lT a ^ ■ \> a - gjD^T^ 



(26) 



where (<f) 1 , F l ) and (V^, A J ', D^) denote the component fields of the chiral 
and vector supermultiplets <3> l and V\ and is the covariant derivative 



D^d^-igjV^a. 



The second Lagrangian of Eq. (124)) . £yang 



-Mills 



(27) 

contains the kinetic terms 



and self-interactions of the vector superfields. It consists in a sum over all 
the vector supermultiplets of the theory, and can be written, for one specific 
vector supermultiplet V, 



C 



Yang-Mills, V 



W a a W a a 



1 



16g 2 L 



W A a W 



(28) 



Similarly to £ cri irab the squared brackets indicate that we only take the co- 
efficient of the corresponding term in the expansion in the Grassmann vari- 
ables. The spinorial superfields W a a and W a a in Eq. ( 128]) . the supersym- 
metric equivalents to the field strength tensor, are related to the quantities 
W a = W a a T a and W & = W A a T a , the latter being given by 



W a = — D ■ D e 29V D a e- 29V and W* = — D ■ D e 
4 4 



2 ° v D A e 29V , (29) 



where D and D denote the superderivatives of Eq. (JSJ) and g the coupling 
constant associated to the gauge group. After expanding W a into component 
fields, one obtains the expression of W a a , 



W a a = -2g 



-iX a „ + 



»7r v t 



CT^CT 



9 n D a - 9 



(30) 



30 



a similar expression existing for W „ a . We recognize in Eq. ( 13 Op the expres- 
sions for the field strength tensors and the covariant derivative in the adjoint 
representation 

IV° = <9 M A° + gf bc a Vft c , 



v* = W - d v v; + gf bc a VX 



(31) 



where fb c a are the structure constants of the gauge group. In the abelian case, 
the expressions for the spinorial superfields can be drastically simplified, 

W a = --D ■ DD a V and W & = -D ■ DD^V , (32) 

and the Lagrangian reduces to 

C um , v = \ [W«W a ] + \ \W a W"} u . (33) 

Finally, the superpotential Lagrangian £ SU perW contains the interactions 
among the chiral superfields derived from the superpotential W(§), a holo- 
morphic function of the chiral superfields. It is given by 



C 



superW 



W($) + W*(&) __ , (34) 



where W^*($t) is the anti-holomorphic function complex conjugate to W($). 
It can be shown that in a renormalizable model the superpotential can at 
most be cubic in the fields, thus taking the form 

W($) = ai$ l + by ¥ $ J + c ijk ¥ & $ fc , (35) 

for some model-dependent parameters a iy b^ and c^. 

As the superpotential is simply a polynomial in the chiral superfields, it 
can be trivially implemented into FeynRules, and the Lagrangian density 
can easily be obtained from the Theta2Component and Thetabar2Component 
functions. The non-trivial part of any implementation of a supersymmetric 
model into FeynRules hence consists in the implementation of £ c hirai and 
^Yang-Miiis- Since these two Lagrangians are however completely fixed by 
gauge symmetry and their form is independent of the actual model under 
consideration, the superspace module of FeynRules comes with some pre- 
defined functions that allow to generate £ c hirai and /^Yang-Mills in an auto- 
mated way. These functions will be described in the next section and are 
summarized in Table El 
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Table [CO Predefined functions related to SUSY Lagrangians. 

CSFKineticTerms [csf ] 

Derives all the kinetic and gauge interaction terms associated to 
the chiral superfield csf. If the function is called without any 
argument, it will sum over the whole chiral content of the theory. 

VSFKineticTerms [vsf ] 

Derives all the kinetic and gauge interaction terms associated to 
the vector superfield vsf. If the function is called without any 
argument, it will sum over the whole gauge content of the theory. 

Superf ieldStrengthL [vsf , alpha, gauge index] 

Calculates the left-handed superfield strength tensor associated 
to the vector superfield vsf. The symbol alpha denotes the free 
spin index, whilst the optional symbol gauge index denotes the 
adjoint gauge index relevant for non-abelian gauge groups. 

Superf ieldStrengthR [vsf , alphadot , gaugeindex] 

Calculates the right-handed superfield strength tensor associated 
to the vector superfield vsf. The symbol alphadot denotes the 
free spin index, whilst the optional symbol gaugeindex denotes 
the adjoint gauge index relevant for non-abelian gauge groups. 

SolveEqMotionD [lag] 

Computes and solves the equations of motion associated to the 
auxiliary D-fields, and then inserts the solution in the Lagrangian 
lag. 

SolveEqMotionF [lag] 

Computes and solves the equations of motion associated to the 
auxiliary F-fields, and then inserts the solution in the Lagrangian 
lag. 
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7.2. Automatic generation of supersymmetric Lagrangians 

The kinetic part of the Lagrangian for a chiral superfield $ can be ob- 
tained automatically in FeynRules from the CSFKineticTerms function. 
As an example, for a chiral superfield implemented as PH30, the Lagrangian 
of Eq. is obtained by issuing the command 

LChiralPhi = CSFKineticTerms [ PHI ] . 

The expression returned by CSFKineticTerms is not automatically expanded 
in terms of the component fields but still expressed in terms of superfields. 
The component-field expression of the Lagrangian can be recovered by apply 
the Theta2Thetabar2Component function to the result, 

Theta2Thetabar2Component [ LChiralPhi ] 

The full Lagrangian £ c hi ra i is obtained by summing over all chiral super- 
fields of the theory. In addition to the function described above which returns 
the kinetic term for a single chiral superfield, FeynRules allows the user to 
obtain directly the complete chiral Lagrangian expressed in terms of super- 
fields via the command CSFKineticTerms. The expression returned by this 
command is equivalent to a sum of terms consisting each in an application 
of the CSFKineticTerms function to a single superfield. The extraction of the 
Lagrangian density can again be achieved via the Theta2Thetabar2Component 
function. Hence, the full Lagrangian £ c hi ra i can be obtained by simply issuing 

Lchiral = Theta2Thetabar2Component [ CSFKineticTerms [ ] ] . 

The supersymmetric equivalents of the field strength tensors can be ob- 
tained automatically in a similar way. The left-handed superfield strength 
tensors W a and W a a associated to a vector superfield V can be called in the 
superspace module of FeynRules via the commands^ 

Superf ieldStrengthL [ V, sp ] , 
Superf ieldStrengthL [V, sp, ga] , 



6 We assume that PHI has been correctly declared in the FeynRules model file. 
7 We assume that the vector superfield V has been declared associated to some gauge 
group in the FeynRules model file. 



33 



in the abelian and non-abelian cases. The symbol sp denotes the undot- 
ted spin index attached to the spinorial superfield while ga is the adjoint 
gauge index relevant for non-abelian gauge groups. Similarly, the abelian 
and non-abelian right-handed superfield strength tensors and Wa a can 
be obtained through 

Superf ieldStrengthR [ V, sp ] , 
Superf ieldStrengthR [ V, sp, ga ] , 

respectively, the only difference with the left-handed case being the variable 
sp which stands this time for a dotted spin index. Note that the spinorial 
superfields defined by Eq. and Eq. (1521) are not hardcoded in the su- 
perspace module of FeynRules, and will be recalculated each time. How- 
ever, the Superf ieldStrengthL and Superf ieldStrengthR functions will 
be evaluated by FeynRules only at the time of the expansion in terms of 
the component fields. 

From the superfield strength tensors we can easily built the kinetic terms 
for vector superfields in an automated way. This is achieved in FeynRules 
by issuing the command 

LV = VSFKineticTerms [ V ] , 

The Super- Yang-Mills Lagrangian of Eq. (128|) can then easily be obtained by 
extracting the 9 ■ 9 and 9 ■ 9 components, 

LSYM = Theta2Component [ LV ] + Thetabar2Component [ LV ] . 

If a model contains several gauge groups, we have to sum over the corre- 
sponding kinetic terms. Similarly to the automatic generation of the kinetic 
terms of chiral superfields, issuing VSFKineticTerms without any argument 
is equivalent to a sum over all possible vector superfields defined in the model. 

At this stage, generating a Lagrangian density for any supersymmetric 
model reduces to an almost trivial task with the help of the functions that 
we have just described. Assuming that a superpotential SP has been defined, 
the full Lagrangian density can be easily implemented into FeynRules as 

Lag = Lchiral + LSYM + LW , 
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where the terms in the sum in the right-hand side are given by 

LC = CSFKineticTerms[ ] , 
Lchiral = Theta2Thetabar2Component [ LC ] , 
LV = VSFKineticTerms[ ] , 
LSYM = Theta2Component [LV] + Thetabar2Component [LV] , 
LW = Theta2Component [SP] +Thetabar2Component [HC [SP] ] . 

The Lagrangian density obtained in this way however still depends on the 
auxiliary F and D fields, which can be eliminated by their equations of 
motion. This can be performed automatically with the help of two functions 
SolveEqMotionD and SolveEqMotionF. Each of them computes, for the in- 
fields and F-fields respectively, the equations of motion directly from the 
Lagrangian, solves them analytically and subsequently inserts the solution 
into the Lagrangian in order to eliminate the auxiliary fields. Using the 
Lagrangian defined previously, the auxiliary fields are eliminated via the 
Mathematica commands 



Lag = SolveEqMotionD [Lag] , 
Lag = SolveEqMotionF [Lag] . 



8. Implementation of the Minimal Supersymmetric Standard Model 

The Minimal Supersymmetric Standard Model (MSSM) is the simplest 
supersymmetric model, resulting from a direct supersymmetrization of the 



Standard Model (SM) (25|,|26|. In this section, we describe its implementation 
in FeynRules in terms of superfields. The corresponding model file can be 
downloaded from the FeynRules website: 



Ihttp : / /f eynrules . phys . ucl . ac . be/wiki/MSSM 



8.1. Gauge groups and representations 

The MSSM is based on the same gauge group as the SM, SU(3) C x 
SU{2)l x U(1)y- The implementation of these three gauge group classes 
slightly differs from the one included into the previous MSSM implemen- 
tation, expressed entirely in terms of the component fields Q. First, we 
associate to each gauge group one vector superfield instead of one gauge 
boson. The abelian factor U(l)y is hence implemented as 
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U1Y == { 

Abelian -> True, 

CouplingConstant -> gp, 
Superfield -> BSF, 

Charge -> Y 

h 

where the vector superfield BSF will be specified below. Secondly, the im- 
plementation of the non-abelian direct factors of the gauge group, SU (2) L 
and SU(3) C , includes a consistent definition of the representation matrices, 
together with the associated index type, related to the representations in 
which one or several chiral superfields of the model liqj. This allows to ex- 
tract the gauge interaction and kinetic terms of the Lagrangian automatically 
with the help of the two functions CSFKineticTerms and VSFKineticTerms. 
For SU(2) L , we only need the fundamental representation, labelled by the 
quantity Ta, and defined together with the associated gauge index SU2D. This 
is implemented as, 

SU2L == { 

Abelian -> False, 

CouplingConstant -> gw, 

Superfield -> WSF, 

StructureConstant -> ep, 

Representations -> {Ta , SU2D} , 

Definitions -> { 

Ta[a__] -> PauliSigma[a]/2, 

ep -> Eps}- 



For the SU(3) C gauge group, both the fundamental and its complex conjugate 
representation are needed, the gauge group being hence defined as, 

SU3C == 
{ 

Abelian -> False, 

CouplingConstant -> gs, 



8 We let FeynRules handle automatically the adjoint representations necessary for the 
vector superfields of the model. We refer to the FeynRules manual for more information. 
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DTerm -> dSUN, 

Superfield -> GSF, 

StructureConstant -> f , 

Representations -> { {T, Colour}, {Tb,Colourb} } 
> , 

where Colour and Colourb are the fundamental and antifundamental rep- 
resentation indices and T and Tb the corresponding representation matrices. 
It is important to note that all antifundamental indices in colour space must 
be replaced by fundamental ones before exporting the Feynman rules to the 
interfaces to Monte Carlo generators, following the conventions of Ref . 0] . 

8.2. Field content 

The Standard Model quarks and leptons are embedded into chiral su- 
permultiplets, together with their squark and slepton partners, which are 
grouped into three generations of six chiral superfields, 



Ql = (3,2,l) , E£ = (3,l,~) , £>*, = (3,1,1) 
L* L = (1,2,-1) , Ek= (1,1,1) , Vk= (1,1,0), 



(36) 



where % stands for a generation index and where we have indicated the rep- 
resentations of the different superfields under the MSSM gauge grourfl The 
component fields included in each superfield can be found in Table [7J For 
completeness, the right-handed neutrino superfield has been introduced, but 
it will be kept sterile, i.e., non-interacting with any other superfield. 

The six chiral superfields of Eq. (1361) are implemented following the in- 



structions given in Section 14.21 while the component fields are implemented 
following the synthax presented in Refs. 0, [22| and the constraints intro- 
duced in Section 14.21 Hence, for each chiral superfield and the associated 
component fields, the hypercharge quantum number and the attached gauge 
indices are specified. This allows the automatic function CSFKineticTerms 
to correctly derive the associated gauge interactions. As an example, the 
weak isospin doublet of quarks Q l L is implemented using the Mathematica 
instructions 

9 It can then be seen that the representation matrices defined in Section [8. II arc enough 
to describe the entire superfield content given in Eq. (|36j) . 
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Tabic 7: The MSSM (s)fcrmion sector, with the representations under the gauge group 
SU(3) C x SU{2)l x U(1)y- The superscript c denotes charge conjugation. 



CSF[1] == { 

ClassName -> QL, 

Chirality -> Left, 

Weyl -> QLw, 

Scalar -> QLs, 

QuantumNumbers -> {Y-> 1/6}, 

Indices -> {Index [SU2D] , Index [GEN] , Index [Colour] }} . 

In our model implementation, we follow a simple naming scheme for the com- 
ponent fields where the names of the Weyl fermionic and scalar components 
are obtained by suffixing w and s to the superfield class name, respectively. 

To preserve the electroweak symmetry from chiral anomalies and in order 
to give masses to both up-type and down-type fermions, the MSSM Higgs 
sector contains two chiral supermultiplets, 

£^ = (1,2,—) , ^ = (1,2,1), (37) 
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Table 8: The MSSM Higgs(ino) sector, with the representations under the gauge group 
SU(3) C x SU(2) L x U(1) Y . 
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Table 9: The MSSM gauge sector, with the representations under the gauge group 
S77(3) c x SU(2) L x U(l) Y - 



each of which consisting in one scalar Higgs SU (2) l doublet and its fermionic 
Higgsino partner, as shown in Table El The superfield Hu couples to up- type 
particles whilst Hd couples to down-type particles. The two chiral superfields 
Hu and Hp are implemented in a similar fashion as presented above. 

The gauge sector is described by three vector superfields associated each 
to one specific direct factor of the gauge group. They lie in the corresponding 
adjoint representation and are singlets under all the other group factors, 



S77(3) C ^$ G = (8,1,0) 
SU(2) L ^$ W = (1,3,0) 
U(1) Y ->> ®b = (1, 1, 0) 



(38) 



These superfields include, in addition to the Standard Model gauge bosons, 
their fermionic partners, the gauginos, as shown in Table [9j Each vector 
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superfield is implemented following the same pattern, with the Indices op- 
tion set to the adjoint index of the relevant gauge group, labelled by SU2W 
and Gluon for SU(2)l and SU(3) C , respectively. As an example, the gluon 
superfield implementation is given by 

VSF[3] == { 

ClassName -> GSF, 

GaugeBoson -> G, 

Gaugino -> gow, 

Indices -> -[Index [Gluon] } 
} • 

8. 3. Lagrangian 

As stated in Section [7J the kinetic and gauge interaction terms of the 
chiral and vector superfields are entirely fixed by gauge invariance. For the 
MSSM, these terms read, 



£sym 
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(39) 



(40) 



4 



D ■ De 29sVa D a e~ 29sVo and W Ga = W%„T a 



In the expressions above, we have introduced the non-abelian vector super- 
fields Vw = $jyfc|crfc, Vq = $G a T a and V G = &G a T a where a k /2 denote the 
generators of the fundamental representation of SU(2)l, a k being the Pauli 
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matrices, and where T a and T a are the generators of the fundamental and 
antifundamental representation of SU(3) C . The gauge coupling constants 
are defined as g', g w and g s and all generation indices are understood. The 
two Lagrangians appearing in Eq. fl39|) are implemented in the model file as 



described in Section IT72| using the automatized functions CSFKineticTerms 
and VSFKineticTerms. 

The interactions among the chiral superfields introduced in Eq. ( 13T)j) and 
Eq. (|37|) are included in the superpotential 



W^mssm = {y^U^Qi-Hu - (y%^ R Q j L -H D - (y%E R L{-Hn 
+ l^Hjj ■ H D , 



(41) 



where y u , y d and y 1 denote the 3x3 Yukawa matrices in flavor space, /i 
the Higgs off-diagonal mass-mixing and the dot products stand for SU (2) in- 
variant products. This superpotential is the most general function satisfying 
renormalizability, the gauge symmetries of the model and .R-parity conser- 
vation. Our conventions regarding the parameters of the model follow the 



SUSY Les Houches Accord (SLHA) [27|, |28| arid the Yukawa matrices must 
then be given flavor- diagonal. However, the fermionic components of the 
superfields given in Eq. f )36|) are gauge-eigenstates and not mass-eigenstates. 
In the model file, we address this issue by implementing the modified super- 
potential 

PUmssm = (y u h U^Qi-Hu - (y d V^ KM ) t3 D l R Qi.H D - {y%^ R V L -H D 
+ IaHjj ■ H D , 

where the hatted Yukawa matrices are flavor-diagonal, following thus the 
SLHA conventions, and where the superfields are gauge-eigenstates. All the 
rotations diagonalizing the gauge-eigenstate basis will be absorbed at the 
component field level rather than at the superfield level, as described in Sec- 
tion |H31 and the role of the CKM matrix appearing in the second term of the 
superpotential is to compensate the only remaining misalignment between 
mass and gauge-eigenstates after these rotations. The Lagrangian associated 
to this superpotential is given by Eq. ( 134)) . 

Finally, the supersymmetry-breaking Lagrangian is obtained by adding 
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explicitly all soft supersymmetry-breaking terms at low-energy, 



C 



Soft 



1 

2 L 



Mi B-B + M 2 W-W + M 3 g-g + h.c. 



m 



.2\i «t /&' 



u HlH u -ml d H\H d (43) 



ij U R ( lL'Hu 



(T 

bH u -H d + h.c. 



The first line of this equation contains the gaugino mass terms, the second 
and third lines the scalar mass terms, m% , 



m?, m?, m| being 3x3 



Q' L 

hermitian matrices in generation space and the fourth and fifth line the bilin- 
ear and trilinear scalar interactions derived from the superpotential. T u , T^, 
and T e are 3x3 matrices in generation space. Similarly to the superpotential 
of Eq. (f42l) . we write the scalar trilinear interactions as, 



t 



3»t d 



H d + h.c. 



(44) 



rather than in their original form, following thus the SLHA conventions. 
Finally, the equation of motions for the auxiliary fields are solved so that 
they are eliminated from the Lagrangian as described in Section 17.21 

8.4- Electroweak symmetry breaking, particle mixings and Dirac fermions 

Due to the strong Yukawa coupling between the superfields Hu, Q\ and 
U R in the superpotential, the electroweak symmetry is radiatively broken to 
electromagnetism and the classical Higgs potential has a non-trivial mini- 
mum. Shifting the neutral scalar Higgs bosons by their vacuum expectation 
values (Vevs), 



and 



v d + K 



(45) 



y/2 V2 ' 

where v u and v d denote the two vevs of the neutral Higgs bosons and h° u and 
h d complex scalar fields, we can extract the mass matrices of the electroweak 
gauge bosons B^ and W*, diagonalize them, and derive the physical mass- 
eigenstates, the photon and the weak bosons and Z^. As in the 
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Standard Model, the transformation rules relating the mass and interaction 
bases are 



where the weak mixing angle 6*^ and the physical masses M z and are 
defined by 

cos 2 ^ = ^^, M Z = -^- and M w = ^ , (47) 
+ # 2 cos 9 W 2 

with = + t^. We define consistently all the mixing parameters in the 
model file, taking the masses as input parameters. The mixing angles and the 
vevs are then dependent on the latter. The rotations are implemented using 
the Definitions option of the particle class. As an example, the SU(2)l 
boson redefinitions are implemented as 

Definitions -> { 

Wi [mu_ , 1] -> (Wbar [mu] +W [mu] ) /Sqrt [2] , 

Wi[mu_,2] -> (Wbar[mu]-W[mu])/(I*Sqrt[2]) , 

Wi [mu_ , 3] -> cw Z [mu] + sw A [mu] 
> , 

where A, Z and W correspond to the model file definitions of the physical 
gauge bosons, containing all the options required in order to have the inter- 
faces to the Monte Carlo codes and to FeynArts working properly (PDG, 
PropagatorType, ...), as presented in Ref. jsf. 

In the Higgs sector, three out of the eight real degrees of freedom of 
the two doublets are the pseudo-Goldstone bosons G ± and G° becoming 
the longitudinal modes of the weak bosons, while the five others mix to the 
physical Higgses, hP, H°, A and H ± . The diagonalization of the scalar, 
pseudoscalar and charged Higgs mass matrices leads to the transformation 
rules 

hi = cos a h° + sin a H° + i cos A + % sin G° , 

h° r ,= - sin a h° + cos a H° + i sin B A - i cos B G° , 
Z , , 48 

H+ = cos/3 H + + sin/3 G + , V ' 

H~[ = sin B H~ — cos 3 G~ , 
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where a is the neutral Higgs mixing angle and the (3 ang le is defined by 



tan/5 = v u /vd- Following the SLHA conventions [27], |28|, all the mixing 
angles are external parameters included in the SLHA block HMIX, and the 
rotations are implemented in the scalar Higgs field class definition, e.g., as 

Definitions -> { 

hus[l] -> Cos [beta] *H + Sin [beta] *GP , 
hus[2] -> (vu + Cos[alp]*hO + Sin[alp]*H0 + 
I*Cos [beta] *A0 + I*Sin[beta] *G0)/Sqrt [2] 

} 

for the H u doublet labelled by hus, where H, AO, hO, HO, GP and GO are the la- 
bels of the physical Higgs and Goldstone bosons. The latter are implemented 
following the synthax presented in Ref. jij and contain again all the options 
required by the interfaces to work properly. 

In the fermionic sector, the mass matrix of the neutral partners of the 
gauge and Higgs bosons is diagonalized through a unitary matrix N which 
relates the four physical (two-component) neutralinos Xi to the interaction- 
eigenstates, 



xl 
xl 
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(iB\ 
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(49) 



Similarly, the mass matrix of the charged partners is diagonalized through 
two unitary matrices U and V relating the interaction-eigenstates to the 
physical (two-component) charginos eigenstates xt according to 



xi 
xi 




and 



Xi 

X2 



u 



iW- 
HT 



(50) 



where the charged winos have undergone the same rotations as the charged 
gauge bosons, 

= ±={Wl =f iWl) . (51) 

The factors of i absorbed in the gaugino definition are required in order to 
obtain real mass matrices. These conventions may seem a bit different from 
the ones specified by the SUSY Les Houches Accord, accounting for factors 
of —i. However, this sign is purely conventional and is related to the sign of 
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the exponentials in the Lagrangian of Eq. (I3"9~j) . The choices of Eq. ( 1491) and 
Eq. (!5Up ensures that the mixing matrices N, U and V are the same as those 
specified in the accord. In the MSSM model file, the real and imaginary 
parts of these mixing matrices are considered as input parameters, following 
the SLHA, while the full matrices themselves are dependent parameters. 
The physical two-component neutralino fields are implemented following the 



synthax of Ref. |22| and the rotations are included in the Definitions option 
of the gauge-eigenstate particle classes. As an example, the wino rotations 
are given by 

Definitions -> { 

wow[s_,l] :> Module [{i}, (Conjugate [UU[i , 1] ] *chmw[s , i] + 

Conjugate [VV [i , 1] ] *chpw[s , i] )/(I*Sqrt [2]) ] , 
wow[s_,2] :> Module [{i}, (Conjugate [UU[i , 1] ] *chmw[s , i] - 

Conjugate[VV[i,l]]*chpw[s,i])/(-Sqrt[2]) ] , 
wow[s_,3] :> Module [{i}, -I*Conjugate [NN [i ,2] ] *neuw[s , i] ] 

> , 

where chmw, chpw and neuw denote the labels of the physical \~ > X + an d X° 
fields, respectively, and NN, UU and VV the mixing matrices. 

The diagonalization of the quark sector requires four unitary matrices, 

tf L ->W £ , d%^U d d%, ui^Vuui, u%^U u u%, (52) 

so that the superpotential of Eq. ( l4Tj) is rotated to a form where the Yukawa 
matrices are diagonal. As a consequence, the charged current weak interac- 
tions become proportional to the CKM matrix 

V C km = V^V d . (53) 

We adopt the standard choice of absorbing these rotations in a redefinition 
of the down-type quark fields alone, 

4, -> VcKudi , (54) 

keeping the up-type quark fields unchanged. This field redefinition is imple- 
mented in FEYN RULES through the Definitions option of the QLw class, 

Definitions -> { 

QLw[s_, 1, ff_, cc_] -> uLw [s ,f f , cc] , 
QLw[s_, 2, ff_, cc_] :> Module [{ff 2} , 
CKM[ff,ff2] dLw[s,ff2,cc] ] 

> , 
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where uLw and dLw denote the (two-component) quark mass eigenstates and 
the CKM matrix is included in the model file through the SLHA blocks VCKM 
and IMVCKM for its real and imaginary parts, respectively, the complete matrix 
being thus an internal parameter. Similarly, the lepton sector is diagonalized 
through the rotations, 

el^V e e* L , eg-H^SS, u* L V„i/« , (55) 

rendering the charged current interactions proportional to the PMNS matrix 

^pmns = VlVv ■ (56) 

We adopt here the choice of absorbing the rotations in a redefinition of the 
neutrino fields alone, 

Vl ~> VpmnsVl , (57) 

leaving the charged lepton fields unchanged^!. The fermion rotations pre- 
sented above legitimate the implementation of the superpotential under the 
form of Eq. f )42|) . It can be checked that the extracted quark and lepton mass 
matrices are indeed diagonal and no further rotation is necessary The PMNS 
matrix is implemented in the model file within the SLHA blocks UPMNS and 
IMUPMNS for its real and imaginary parts, respectively, while again, the entire 
matrix is considered as an internal parameter. On a similar fashion as for 
the quark sector, the implementation of the field redefinitions is given by 

Definitions -> { 

LLw[s_, 1, ff_] :> Module [{ff 2}, 
PMNS [f f , f f 2] *vLw [s , f f 2] ] , 

LLw [s_ , 2 , f f _] ->eLw [s , f f ] 
> , 

where LLw is the gauge-eigenstate field and vLw and eLw the mass-eigenstate 
ones, the right-handed components being unchanged. 

In the scalar sector, we define the super-CKM and super-PMNS bases 
as the bases in which the scalar fields undergo the same rotations as their 
fermionic counterparts. However, the fermion and sfermion fields can be 



10 There are many different modefs accounting for neutrino masses. Since there is thus 
not any unique way to define neutrino mass terms, we are not including them in the imple- 
mented Lagrangian and we only account for neutrino mixings through field redefinitions. 
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misaligned due to possible off-diagonal mass terms in the supersymmetry- 
breaking Lagrangian £soft ; and four additional rotations R u , R d , R e and R v 
are (in general) required, 
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(58) 



where the flavor-eigenstates are denoted explicitly and where the physical 
eigenstates are mass-ordered, from the lightest to the heaviest. Again, the 
four rotation matrices are implemented after splitting their real and imag- 
inary parts following the SLHA-2 conventions. The physical scalar fields 
are implemented following the synthax of Ref. ji[ and contains all the in- 
formation needed by the various FeynRules interfaces, including the field 
redefinitions in the Definitions option. As an example, the redefinition of 
the scalar component of the superfield U l R is given by 

Definitions -> { 

URs[ff_, cc_] :> Module [{ff 2}, subar [f f 2 , cc] *RuR[f f 2,f f ] ] 
> , 

where RuR refers to the three last columns of the mixing matrix R u and su 
denote up-type squark mass-eigenstates. 

Finally, we define Dirac representations ip f° r the fermions, because most 
generators work in terms of Dirac fermion rather than Weyl fermions, the 
former being required at the Monte Carlo generator level. In terms of the 
Weyl fermions introduced above, the Dirac fermions read 
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The four-component Dirac and Majorana fermions are implemented following 
the synthax of Ref. [3], where we associate a four-component fermion to its 
Weyl components through the WeylComponents option of the particle class 



221 ] . As an example, the charged lepton is implemented as 



F[l] == { 

ClassName->l , 
Self Conj ugate->False , 
Indices->{Index [GEN] } , 
Flavor Index->GEN, 
WeylComponents->-[eLw , ERwbar} , 



where the dots stand for additional options such as those required by the 
Monte Carlo tools. 

The FeynRules function WeylToDirac allows to perform the expansion 
of the Weyl fermions in terms of the Dirac fields. However, we first need to 
address the issue related to the antifundamental color representation which 
the right-handed quark fields lie in. This is mandatory in order to have one 
single color representation for the Dirac fields, i.e., the fundamental one, 
and for the interfaces to the various tools linked to FeynRules to work 
properly. Denoting T and T the fundamental and antifundamental color 
representations, and using the property T = — T*, the problem is solved by 
implementing the instructions, 

Colourb = Colour , 

Lag = Lag /. { Tb [a_ , i_ , j_] ->-T [a, j , i] } , 

where Lag denotes the MSSM Lagrangian. Then, since Dirac fermions are 
defined only for physical particles, we start with an expansion of the SU(2) l 
multiplets in terms of their components, 

Lag = Expandlndices [ Lag , FlavorExpand -> {SU2W, SU2D} ] , 

before eliminating the Weyl fermions from the Lagrangian, 

Lag = WeylToDirac [ Lag ] . 
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The obtained Lagrangian is now suitable for the calculation of the Feyn- 
man rules through the function FeynmanRules or to be exported to Fey- 
nArts with the function WriteFeynArtsOutput or to any Monte Carlo tools 
linked to FeynRules via the functions WriteCHOutput, WriteMGOutput, 



WriteSHOutput, WriteUFG [29| or WriteWOOutput (5| 



8.5. Validation 

Our implementation has been validated against the (public) MSSM im- 
plementation in the current FEYN RULES version 1.4.x. The latter is based on 
component fields and Dirac fermions and has been validated both against the 
literature and the built-in MSSM implementations in various Monte Carlo 
tools (I)]. We have checked analytically that the Feynman rules obtained with 
FeynRules using the superfield MSSM implementation correspond to those 
found in the literature and verified that 

FeynmanRules [ LagSF - LagComponents ] 

provides an empty list of Feynman rules, ensuring that the two Lagrangians 
are equal. In the expression above, LagSF denotes the Lagrangian obtained 
using the superspace module of FeynRules and LagComponents is the one 
available in the public version of the model. 

8.6. Gauge choice 

In Eq. (I48p . we have introduced the Goldstone bosons absorbed by the 
weak gauge bosons to get their longitudinal polarization. Although absent in 
calculations performed in unitarity gauge, they must be included if another 
gauge is used. In addition, a gauge-fixing Lagrangian must be included in the 
model, as well as ghost fields associated to the gauge bosons, together with 
the corresponding Lagrangian constructed from the BRST formalism. These 
two Lagrangians are included in our MSSM implementation in the specific 
choice of the Feynman gauge. 

The gauge-fixing Lagrangian is most conveniently derived after rewriting 
the Higgs doublets in terms of eight real scalar fields 0* with i = 1, ... ,4. 
The Higgs doublets become 





(<j>l + ic 


*2 




{<f>l + ic 












K<t>d + H 





Hd= 7i (i + ti) - (0, °' °' °' ^ ^ 



(60) 
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and the gauge-fixing Lagrangian is given by 



£gf — — r (g b G b + G^G^ + G^G^ 



(61) 



with 



/~rfe 



W 





- 9 



l\-\(H u ))-H u -g'^-(H d ))-H d , 
9v,y - l(<T k Hu)) ■ H u - g w (^((T k H d )j ■ H d , 



(62) 




,u 



where the dot product stands for the scalar product of the field space intro- 
duced above. Comparing with Eq. ( |4"5]) . we derive the only non-zero vevs as 
(4>l) = v u and (4> l d ) = v d . 

The Fadeev-Popov ghost Lagrangian is derived from the gauge variation 
of the gauge- fixing functions and reads, for the S-ghost ub, the H^-ghost uw 
and the gluon ghost u g , 

£ g host = -ubO^ub - u w d^D^u w - Ugd^D^Ug + C ghosttlf> , (63) 

where the adjoint representation indices are understood and the covariant 
derivatives are taken in the adjoint representation, 



D^ w = d^ w + gJ jk Wlu k w , DXg = dX g + 9sf\c9lu c g . (64) 



with ua = ub, u\y, u^r, u^, the corresponding real representation matrices 
being T A = —iY, — ia 1 /2, —ia 2 /2, —ia 3 /2 and the dot product stands for 
the scalar product defined in the basis {4> z u , (fi d } introduced above. The ghosts 
related to the physical gauge bosons are obtained by rotating those related 
to the gauge-eigenstates parallel to the latter. 

The MSSM implementation in FeynRules contains the two Lagrangians 
of Eq. flSTj) and Eq. fl63|) and is hence fully expressed in Feynman gauge, the 
unitarity gauge being recovered by removing all the ghosts and Golstone 
bosons from the Lagrangian. 



Finally, the scalar piece of the ghost Lagrangian is given by 



Aw = u A [(T A (H u )) ■ {T A '[(H U ) + H u }) + 

{T A (H d )) ■ {T A '[(H d ) + H d ])]u A > , 



(65) 
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9. Conclusions 

In this paper we presented a superspace module for the FeynRules 
package, and the whole module is distributed together with the the Feyn- 
Rules package starting from version 1.6.x. The module allows to perform 
computations of superspace quantities involving Grassmann variables and is 
well suited for the implementation of Lagrangians in terms of superfields. 
The package hence allows to implement supersymmetric models into Feyn- 
Rules directly in terms of superfields, the superfield expressions can then 
be expanded into component fields and the corresponding coefficients of the 
Grassmann variables extracted in an automated way. Furthermore, since the 
only piece of a supersymmetric Lagrangian that is not fixed by supersym- 
metry and gauge invariance is the superpotential, the package allows for an 
automatic generation of all the kinetic terms and gauge interaction terms 
for the superfields, thus reducing the implementation of a supersymmetric 
Lagrangian to the almost trivial task of writing down the superpotential. 
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